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ABSTRACT 

The early universe hosted a large population of low-mass virialized "minihalos," that were not 
massive enough to form stars on their own. While most minihalos were photoevaporated by ionizing 
photons from star- forming galaxies, these galaxies also drove large outflows, which in some cases would 
have reached the minihalos in advance of ionization fronts. In the previous papers in this series, we car- 
ried out high-resolution, three-dimensional adaptive mesh refinement simulations of outflow-minihalo 
interactions that included non-equilibrium chemistry, radiative cooling, and turbulent mixing. We 
found that, for a fiducial set of parameters, minihalos were transformed into dense, chemically ho- 
mogenous stellar clusters. Here we conduct a suite of simulations that follow these interactions over 
a wide range of parameters including minihalo mass, minihalo formation redshift, outflow energy, 
outflow redshift, distance, concentration, and spin. In almost all cases, the shocked minihalos form 
molecules through nonequillibrium reactions and then cool rapidly to become compact, chemically- 
homogenous stellar clusters. Furthermore, we show that the unique properties of these clusters make 
them a prime target for direct study with the next generation of telescopes, and that there are many 
reasons to suspect that their low-redshift counterparts are the observed population of halo globular 
clusters. 

Subject headings: galaxies: formation - galaxies: high-redshift - star clusters: general - globular clus- 
ters: general - shock waves - galaxies: abundances 



1. INTRODUCTION 

In the modern theory of cosmic evolution, structures 
form hierarchically through accretion and mergers {e.g. 
White & Rees 1978; White & Frenk 1991; Kauffmann et 
al 1993; Cole et al 2000; Bower et al 2006). A generic 
feature of this history, then, is the presence of a high- 
redshift population of gravitationally-bound structures 
that were much less massive than present-day galaxies. 
At virial temperatures below 10 4 K, atomic transitions of 
hydrogen and helium are not excited, and the gas must 
radiate energy via molecular transition lines or dust emis- 
sion. Although some H2 left over from recombination 
was able to cool the earliest structures {e.g. Abel et al 
2002; Bromm et al 2002; Stacy et al 2010; Turk et al 
2009), the resulting 11.2-13.6 eV background emission 
from these objects {e.g. Haiman et al 1997,2000; Ciardi 
et al 2000; Sokasian et al 2004; O'Shea & Norman 2007) 
quickly dissociated the already low abundance of these 
primordial molecules (Galli & Palla 1998). In fact, even 
if some small fraction of molecules survived, it is unlikely 
to have had a large effect on the structure of these low- 
mass objects {e.g. Whalen et al 2008a; Ahn et al 2009). 
The result of this inefficient cooling then, was a large 
population of primordial, atomic "minihalos" that were 
unable to form stars until triggered by some exterior in- 
fluence. 

A prime candidate for inducing star formation in these 
objects is galaxy outflows, powered by core-collapse su- 
pernovae and winds from massive stars. Such outflows 
have been found originating from a variety of galaxies, 
from dwarfs to massive starbursts, over a wide range 
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of redshifts {e.g. Lehnert & Heckman 1996; Franx et al 
1997; Pettini et al 1998; Martin 1999; Heckman et al 
2000; Veilleux et al 2005; Rupke et al 2005; Chung et al 
2011). Furthermore, theoretical studies have argued that 
many of these early galaxies represent the tail end of a 
large population of low-mass starbursts that occurred be- 
fore reionization (Scannapieco, Ferrera, & Madau 2002; 
Thacker, Scannapieco, & Davis 2002). While ionizing 
photons generated by these early galaxies would have 
lead to the photoionization of large regions of space, these 
photons are also easily trapped behind outflows as they 
sweep up intergalactic gas (Fujita et al 2003). This sug- 
gests that many regions in the intergalactic medium may 
have be impacted by outflows before they were ionized. 

In the first two papers in this series (Gray & Scanna- 
pieco 2010, hereafter Paper I; Gray & Scannapieco 2011, 
hereafter Paper II) we carried out high-resolution, three- 
dimensional adaptive mesh refinement simulations of the 
interaction between a typical primordial minihalo and 
high redshift galaxy outflow, capturing all the important 
physical processes in detail. In Paper I we looked at 
the influence of a 14-species primordial non-equilibrium 
chemical network with associated cooling terms, and 
showed that much of the minihalo baryonic matter is 
removed from the dark matter halo and compressed into 
several dense clumps embedded in a ribbon of gas. The 
inclusion of a dissociating UV background, while altering 
the final molecular abundances, had no effect on the final 
distribution. In Paper II we extended this study by in- 
cluding the effects of metal-line cooling and turbulence. 
Turbulence had the primary effect of mixing metals into 
the primordial cloud from the outflow, which enriched 
the cloud to Z « 10 -2 Z . While metals provided an- 
other source of cooling, we found that they did not have 
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a significant impact on the final structure of the cloud, 
which was very similar to the that found in Paper I. In 
both studies, dense ~ 10 5 M clusters were formed. 

Additionally, these clusters had many of the properties 
of present day halo globular clusters. In particular, the 
halo globular cluster mass function is well defined by a 
Gaussian with a mean mass of 10 5 M and dispersion of 
0.5 dex (e.g. Armandroff 1989), which spans the range 
of clusters formed in our studies. Furthermore, while the 
stars in a typical globular cluster have a small scatter 
in their metallicities, the metallicity distribution func- 
tion of halo globular clusters is also well defined by a 
Gaussian with a mean value of [^] ~ -1.6 and a dis- 
persion of 0.3 dex (Zinn 1985; Ashman & Bird 1993). 
The clusters we formed in Paper II were enriched homo- 
geneously via turbulent processes to approximately this 
mean value. Finally, our models naturally reproduced 
observations that suggest that globular clusters do not 
reside within dark matter potentials (Moore 1996; Con- 
roy et al 2010), as the outflows always removed the halo 
gas from their associated dark matter. 

These many tantalizing clues hint that outflow- 
minihalo interactions may in fact have provided the un- 
known mechanism that formed the observed population 
of halo globular clusters. However, a proper assessment 
of this connection is not possible without understanding 
how robust the mechanism is to the many likely vari- 
ations in the properties of high-redshift minihalos and 
galaxy outflows. Thus, here we turn our attention to 
the range of interactions that lead to the formation of 
compact stellar clusters, carrying out a detailed param- 
eter study in which we vary a wide range of properties, 
including minihalo mass, minihalo formation redshift, 
outflow energy, outflow redshift, distance, concentration, 
and spin. With these results we are able to asses in 
detail the viability of this mechanism in producing halo 
globular clusters, as well as make a series of high-redshift 
predictions to be verified with the next generation of tele- 
scopes. 

The structure of the paper is as follows. In Section §2 
we describe the model setup and parameters, and §3 we 
present the results of these simulations. Finally, in §4 we 
present simulated images of narrow band Lya observa- 
tions from the next generation of ground based telescopes 
and wide band imaging with the James Webb Space Tele- 
scope (JWST) and relate our models to present-day halo 
globular clusters. Our conclusions are given in §5. 

2. MODEL FRAMEWORK 
2.1. Numerical Methods 

All simulations were performed with FLASH ver- 
sion 3.1.1, a multidimensional adaptive mesh refinement 
(AMR) hydrodynamics code (Fryxell et al 2000) which 
solves the Riemann problem on a cartesian grid using 
a directionally-split Piecewise-Parabolic Method (PPM) 
solver (Colella & Woodward 1984; Colella & Glaz 1985; 
Fryxell, Miiller, & Arnett 1989). In addition to mak- 
ing use of the efficient multigrid self-gravity package in- 
cluded with FLASH (Ricker 2008), we have developed a 
variety of other physics packages necessary for simulating 
outflow- minihalo interactions. 

First, we have implemented a version of the non- 
equilibrium chemistry network from Glover & Abel 



(2008). This tracks the full set of two-body reactions that 
determine the evolution of three atomic hydrogen species 
(H, H + , and H~), three atomic helium species (He, He + , 
and He ++ ), three atomic deuterium species (D, D + , and 
D~), two states of molecular hydrogen (H2 and H^), two 
states of hydrogen deuteride (HD and HD + ), and elec- 
trons (e~). However, we have neglected the contribu- 
tion from three-body reactions, since they only become 
important at densities of n > 10 8 cm -3 (e.g. Palla et 
al 1983), and neglected the contribution from reactions 
concerning molecular deuterium (D2), because the cool- 
ing from D2 and is negligible (Glover et al 2008). 
While we included UV dissociation rates from Glover & 
Savin (2009) to capture the effect of a dissociating back- 
ground, we discovered that for all reasonable choices of 
background intensity, they did not effect the evolution of 
these interactions significantly (Paper I). Thus we did not 
include a dissociating background in the present study. 
Further details of this implementation and the impact of 
the UV background are given in Paper I. 

Secondly, we have implemented molecular and atomic 
cooling rates. At temperatures below 10 4 K cooling is 
dominated by both molecular line transitions from H2 
and HD and metal line transitions. We have imple- 
mented the cooling rates from Glover & Abel (2008) for 
the collisional excitation between H2 and a variety of 
atomic species as well as the Lipovka et al (2005) rate 
between HD and H. Above 10 4 K, primordial cooling is 
dominated by atomic hydrogen and helium, whose cool- 
ing rates are calculated using CLOUDY (Ferland et al 
1998) assuming Case B recombination. The metal-line 
cooling function is taken from Wiersma et al (2009) via 
a table lookup, assuming solar abundances ratios, scaled 
by the local metal abundance. 

Finally, we have implemented a buoyancy and shear 
driven model of turbulence that extends the two equa- 
tion K — L model of Dimonte & Tipton (2006; see also 
Chiravalle 2006), where K represents the turbulent ki- 
netic energy and L is the eddy length scale. This model 
reproduces the effect of three primary instabilities: the 
Ray leigh- Taylor instability that arises when a low den- 
sity fluid supports a high density fluid under an accel- 
eration, the Richtmyer-Meshkov instability, which ap- 
pears when a shock interacts with a fluid with a different 
impedance (such as a density gradient), and the Kelvin- 
Helmholtz instability, which occurs between two fluids 
that shear in a direction perpendicular to their inter- 
face. This model also includes diffusion terms for both 
the molecular species and metal abundances, which al- 
lows for metal mixing. Further details of this package 
are given in Paper II. 

As in Paper I and Paper II, we assume a standard 
ACDM cosmology with h = 0.7, Q = 0.3, Q A = 0.7, 
and = 0.045 (e.g. Spergel et al 2007), where h is the 
Hubble constant with units of 100 km s _1 Mpc -1 , Qq, 
^a, are the total matter, vacuum, and baryonic mat- 
ter densities, in units of the critical density. With our 
choice of h, the critical density is p cr i t = 9.2 x 10 -30 g 
cm -3 . 

2.2. Minihalo 

In all runs the minihalo was comprised of a dark matter 
halo and a metal- free, neutral atomic cloud made up of 
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76% hydrogen and 24% helium by mass. The total mass 
of the minihalo was defined as M c = M6 x 10 6 M , and 
we assumed that the dark matter and gas have collapsed 
by a redshift z c at which time the object had a mean 
overdensity of A = 178 (e.g., Eke et al 1998). This leads 
to a mean density of p c = A (1 + ^c) 3 Pcrit in the 
cloud. With these choices the virial radius of the cloud 
is 



R r = 0.3M 



and its virial velocity is 
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The radial profile of the minihalo was taken from 
Navarro et al (1997): 



p(R) = 



Q p c 



cx(l + cx) 2 3F(c) 



g cm 
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where c is the halo concentration parameter, x = R/R c , 
and F(t) = ln(l + 1) — As the minihalo collapses, 
we assume that the gas is isothermal and heated to its 
virial temperature of T c = 720M 6 2/3 [(l + z c )/10] K, which 
produces an isothermal gas density distribution in the 
CDM potential well of 
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where v 2 sc (R = xR c ) = 2v 2 [F(cx)+cx/ (1+cx)] [xF(c)]- 1 
is the escape velocity of a particle at a distance R from 
the halo center, v 2 sc (0) = 2v 2 c/F(c), and c is the halo 
concentration factor which has a fiducial value of 4.8 
(Madau et al 2001). The central gas density of the cloud 
is then 
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where A = 2c/F(c) and t = cx. Outside the virial radius 
of the cloud, we assume that the gas is at the same tem- 
perature as the minihalo and in hydrostatic balance with 
it. This leads to a density distribution of 
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p(R > R c ) = p(R c )e R R c gm cm 
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where R is the radius, R c is the virial radius, and Ro = 
GM c m p /kt)T, where M c and T are the mass and virial 
temperature of the cloud. 

We begin each simulation with the halo in hydrostatic 
equilibrium using a two part gravity scheme to account 
for the self-gravity of the gas and the dark matter halo. 
The self-gravity was handled by the efficient multigrid 
Poisson solver (Ricker 2008) which is included in FLASH. 
The dark matter term was handled by first calculat- 
ing the total gravitational acceleration due to the to- 
tal (dark matter + gas) mass distribution (via Eqn. 2]) 
and subtracting off the term for the initial gas configura- 
tion before adding the results of the self-gravity calcula- 
tion. Thus, initially when the cloud is pressure supported 
against collapse the gas-only and self-gravity terms can- 
cel. Outside the acceleration due to the dark matter was 
calculated as —GM C /R 2 cm s -2 , and again the gas was 



taken to be in initial hydrostatic equilibrium according 
to Eq. H 

2.3. Outflow 

The galactic outflow was modeled as a Sedov- Taylor 
blast wave. We assumed that the minihalo sits at a 
(physical) distance R s and that the shock moves with 
a velocity of 
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where S44 = 5/44, and S is defined as the ratio of the 
density of the gas compared to the mean density at that 
redshift, eE^ is the input energy of the shock wave where 
e is the wind efficiency and £"55 is the total SNe energy in 
units of 10 55 ergs, and z s is the redshift when the outflow 
reaches the halo (Scannapieco et al 2004). The postshock 
temperature is T s = 1.4 x 10 5 (v s /100 km s -1 ) 2 K. By 
the time the outflow reaches R s it has entrained a mass 
of 
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and has a surface density of 
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The outflow is assumed to consist of completely ion- 
ized hydrogen and helium with the same relative abun- 
dances as the primordial gas, but enriched with metals 
from the supernovae that drove the material out of the 
host galaxy. To determine the metal abundance, we fol- 
lowed the estimate from Paper II and Scannapieco et 
al (2004) and assumed that each supernova, whether 
a core-collapse or more-exotic pair-instability supernova 
(Woosley & Weaver 1995; Heger & Woosley 2002), gen- 
erates 2 M of metals per 10 51 ergs of energy. If half 
of these metals are funneled into the outflow, we can ex- 
pect a total mass of metals of M meta i = 10 4 £^55 M . We 
therefore initialized the shock with an initial abundance 
of 



Mmetai 



M, 



metal 



10 4 E 55 M Q 



M 



s,total 1.4 X 10 6 (5 44 (^t) R 3 S 



(10) 



The lifetime of the shock can be estimated from <j s = 
Ppost^postts, where a s is the surface density of the out- 
flow, Ppost is the post shock density, v post is the post 
shock velocity, and t s is the shock lifetime. After a pe- 
riod tfuih which is nominally defined as 1.5 Myrs, the 
shock was allowed to taper off by slowly lowering the 
density and raising the temperature so that the pressure 
stayed constant. This prevents the courant time step 
from becoming exceedingly small long after the shock 
had passed over the minihalo. For more setup details see 
Paper I and Paper II. 

3. PARAMETER STUDY 

Table[T] summarizes the range of parameters studied in 
our simulations. Each case is given a name based on the 
value of the parameter that was changed from the fidu- 
cial case. OFID is the fiducial model from Paper II, in 
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TABLE 1 

Summary of Study Parameters. M6 is the minihalo mass in units of M6 = M c /10 6 Mq, E55 is the energy of the shock in 

UNITS OF E/10 55 ERGS, R s IS THE DISTANCE BETWEEN THE GALAXY AND THE MINIHALO IN UNITS OF (PHYSICAL) KPC, Z c IS THE RED SHIFT 
AT WHICH THE HALO VIRIALIZES, Z s IS THE REDSHIFT AT WHICH THE SHOCK REACHES THE MINIHALO, A' IS THE SPIN PARAMETER (SEE £|3. 6|) , 
C IS THE CONCENTRATION PARAMETER, Z IS THE METAL ABUNDANCE OF THE OUTFLOW IN UNITS OF SOLAR METALLICITY (Zq) , AND RES IS 

THE MINIMUM RESOLUTION OF EACH SIMULATION IN UNITS OF PC. 



Name 


M 6 


E55 


R s 


Z c 


z s 


A 7 


c 


Z 


Res 


Notes 








(kpc) 










(z©) 


(PC) 




OFID 


3 


10 


3.6 


10 


8 





4.8 


0.12 


9.14 


Original Fiducial 


PM10/NFID 


10 


10 


3.6 


10 


8 





4.8 


0.12 


19.85 


New Fiducial 


PM03 


0.3 


10 


3.6 


10 


8 





4.8 


0.12 


5.06 




PM30 


30 


10 


3.6 


10 


8 





4.8 


0.12 


19.85 




PE1 


10 


1 


3.6 


10 


8 





4.8 


0.06 


19.85 




PE5 


10 


5 


3.6 


10 


8 





4.8 


0.06 


19.85 




PE20 


10 


20 


3.6 


10 


8 





4.8 


0.24 


19.85 




PE30 


10 


30 


3.6 


10 


8 





4.8 


0.36 


19.85 




PZC8 


10 


10 


3.6 


8 


8 





4.8 


0.12 


19.85 




PZC15 


10 


10 


3.6 


15 


8 





4.8 


0.12 


19.85 




PZS10 


10 


10 


3.6 


10 


10 





4.8 


0.12 


19.85 




PR21 


10 


10 


2.1 


10 


8 





4.8 


0.60 


19.85 




PR66 


10 


10 


6.6 


10 


8 





4.8 


0.019 


19.85 




PR120 


10 


10 


12 


10 


8 





4.8 


0.003 


19.85 




PSPZ 


10 


10 


3.6 


10 


8 


0.023 


4.8 


0.12 


19.85 


Angle = 0° 


PSPN 


10 


10 


3.6 


10 


8 


0.023 


4.8 


0.12 


19.85 


Angle = 90° 


PC32 


10 


10 


3.6 


10 


8 





3.2 


0.120 


19.85 




PC73 


10 


10 


3.6 


10 


8 





7.3 


0.12 


19.85 





which it was labeled MRWT. In all simulations the wind 
efficiency was e = 0.3, and the mean over density of the 
medium between the outflow and the minihalo was 44 
(J44 = 1). Each simulation was performed in a rectangu- 
lar box with the x-axis twice the size of the y- and z-axes. 
The minihalo was centered at [0,0,0] and the shock origi- 
nated from the left x boundary. The base grid was taken 
to have 16 by 8 by 8 cells in the x, and z directions, 
and l max = 5, such that up to 4 additional levels of refine- 
ment were added in regions with significant pressure or 
density structure. In Paper I and Paper II we varied the 
maximum resolution in our simulations, and showed that 
such medium-resolution (256 x 128 x 128 effective) sim- 
ulations were able to faithfully reproduce the outcome 
of shock-outflow interactions. In the interest of speed 
we also included a forced-derefinement routine that re- 
quired regions to derefine if their density was less than 
3.0xl0 -26 g cm -3 . This was implemented outside 393 
pc from the center of the halo and after 7 Myrs from the 
beginning of the simulations, such that it only impacted 
low-density regions far from the evolving cloud. 

Each simulation was run to a time at which the shock 
completely overran the minihalo, which was typically on 
the order of a few Myrs. However, the clouds are ex- 
pected to evolve over hundreds of Myrs before reaching 
a final configuration, a time that is much longer than we 
were able to run our FLASH simulations. To get around 
this limitation, we transformed the final mass distribu- 
tion in our 3-D simulations into a 1-D ballistic problem, 
as described in Paper I. 

Here, the mass distribution at the final output from 
each simulation was divided into 100 evenly spaced bins 
along the x-axis, and the mass of each bin was calculated 
by summing together the gas within a cylinder with a ra- 
dius of 10 21 cm and a length set by the bin spacing. Next 
each bin was converted into a particle, whose initial po- 
sition was located at the center of mass of the summed 



density distribution, and whose initial velocity was cal- 
culated from conservation of momentum. These particles 
were then evolved ballistically using a leap-frog method. 
At every time step, the acceleration of each particle was 
calculated from the self-gravity from all other particles 
as well as from the gravity from the dark matter halo. 
If a particle moved past the one in front of it, they were 
merged by adding together their masses and calculating 
a new velocity from momentum conservation. The out- 
comes from the different runs were then compared with 
each other at a time of 200 Myrs after the end of each 
simulation. 

3.1. Effect of Halo Mass 

The outflow-minihalo interaction goes through several 
evolutionary stages, which are illustrated in Figure [TJ As 
the outflow impacts the front of the minihalo, the atomic 
gas is ionized but then cools rapidly. As the temperature 
approaches w 10 4 K, molecules begin to form, catalyzed 
by H~ and H + . Since the shock moves faster through the 
less dense regions of the minihalo, molecule formation 
predominately occurs in the gas surrounding the center 
of the impacted minihalo (top left panel). When the 
outflow meets at the back of the minihalo, however, a 
second shock is driven backwards toward the center of 
the collapsing cloud. This promotes molecule formation 
in the center of the minihalo, which continues to collapse 
as it further loses pressure support (top right panel). As 
the cloud becomes more elongated (bottom left panel), 
the shear layer between the outflow and cloud becomes 
more and more turbulent, mixing metals into the cloud. 
The final state of the simulation (bottom right panel) is 
similar to that seen in the runs in Paper I and Paper 
II: the gas from the minihalo has been moved out from 
the dark matter halo and formed a dense, cold ribbon of 
material stretched along the x-axis. 

As the total mass in the impinging wind is relatively 



Fig. 1. — Rendered snapshots of the logarithm of the H2 density for run PM10/NFID. Colors show logarithmic contours of H2 density 
from 1CT 32 to 1CT 26 g cm . Top Left: t = 12.2 Myrs, showing the outflow completely surrounding the minihalo and H2 formation starting 
at the front of the cloud. Top Right: t = 15.7 Myrs, showing the minihalo as it is collapsing and during which H2 is increasing towards the 
center of the cloud. Bottom Left: t = 19.2 Myrs, showing both the cloud as it is stretched and the H2 that is then found throughout the 
dense ribbon. Bottom Right: t = 22.0 Myrs, showing the final state of run PM10 and the uniformly cold ribbon of gas. 



small, the minihalo is the primary source of gas for the 
final dense stellar clusters. Figure [2] shows the impact 
of changing this mass, illustrating the final state of runs 
with minihalo masses ranging between 0.3 to 30 xlO 6 
M©. In all cases the gas from the minihalo forms dense 
cold clumps embedded in a ribbon of diffuse gas, stretch- 
ing away from the dark matter halo. The larger the halo 
mass, the longer the ribbon. In the larger halos, a very 
slight metallicity gradient is observed as some of the gas 
in the center of the halo is not as enriched as the sur- 
rounding gas. However, every part of the ribbon is en- 
riched to well above 1O _3 Z , and the majority of the gas 
is enriched to a nearly constant value of « 10~ 2 Zq. 

Figure [3] shows the evolved distribution of stellar clus- 
ters. The relative sizes of each cluster corresponds to 
their final mass, with larger symbols representing larger 
clusters. As expected, the more massive the initial halo, 
the more massive the final clusters. In the case of the 
largest mass minihalo, the outflow is not strong enough 
to remove all the gas. Instead it leaves behind a rather 
sizable cluster at the center of the dark matter halo, al- 
though even more massive clusters are formed behind it. 
In all other cases, the outflow effectively removes all of 
the gas and forms several dense clusters. Since the final 
cluster sizes are larger in the PM10 model relative to the 
Mq = 3 run, but 1O T M0 solar mass minihalos are nev- 
ertheless quite common at high-redshift, we have used it 
as our new fiducial model (NFID) and maintained this 
value in the simulations discussed below. 



3.2. Effect of Shock Energy 

The shock plays three important roles in the evolution 
of these interactions. First, it accelerates the minihalo 
gas with enough momentum to move it out of the dark 
matter halo, second, it provides enough energy to ionize 
the gas, trigging the non-equilibrium chemical reactions 
that provide an avenue for H2 and HD formation and 
cooling, and finally, the shock brings in the metals that 
enrich the post-shock star-forming medium. The single 
most important parameter in determining each of these 
effects is the energy driving the shock, £^5, as this not 
only sets the shock velocity (Eq. [7]), but determines the 
total associated metal mass and abundance (Eq. [TO]) . 

The fiducial estimate for E55 is taken from Scannapieco 
et al (2004) where they estimate the energy in the outflow 
of a high-redshift starburst galaxy. First, consider a a 
young galaxy of total mass M g = 10 9 M that has a 
virial temperature greater than 10 4 K and is therefore 
allowed to cool via atomic hydrogen. If 10% of the gas is 
converted into stars and 10 51 ergs of energy is produced 
by SNe for every 30 M of stellar material for massive 
stars (Pop III, Heger & Woosley 2002) then the expected 
total energy given off is £^5 ~ 50Mg, where Mg is the 
mass of the galaxy in units of 10 9 Mq. In simulations 
by Mori et al (2002) the wind efficiency e is found to 
be ~ 0.3 for a 2.0 x 1O 8 M starbursting galaxy. Using 
these numbers, the fiducial value for E55 is obtained. To 
get an idea of how this parameter alters the evolution 
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Fig. 2. — Comparison of final outputs from runs with varying halo masses. The first column shows logarithmic density contours from 
1CT 26 to 1CT 21 g cm 3 , the second column shows logarithmic metallicity contours between 10 4 and 10 3 Zq, and the third column 
shows logarithmic temperature contours between 10 and 10 8 " 6 K. The first row shows results from run PM03, the second row shows results 
from OFID, the third row shows results from run PM10, and the last row shows results from run PM30. In all cases the outcome is similar. 
The minihalo gas is moved out from the dark matter halo, stretched along the a>axis, and enriched homogeneously. 



of the minihalo, we allow this value to range between 
1 < < 30. 

Figure 2] shows the results of altering this key param- 
eter. Each row corresponds to a different value of E^, 
which increases from top to bottom as 1, 5, 10, 20, and 
30, respectively. As the shock energy increases, the mini- 
halo evolution changes dramatically. Instead of stretch- 
ing the gas into a ribbon as in the fiducial case, £^5 = 10 
case, the stronger shocks crush the cloud into a single 
small, dense cluster. Interestingly, even though these 
high-energy interactions happen quickly, there is still 
enough time to mix in the metals from the shock. In all 
cases the metallicity is ~ 10 _2 Zq, a value that is roughly 
constant across models in part because the higher energy 
shock models have higher initial metal abundances. On 
the other hand, the lower the shock energy, the more the 
cloud is stretched into a diffuse ribbon in which smaller 
clusters are embedded. 

Over time, much of the surrounding gas merges with 
the initial cluster, which is similar to the evolution we 
found in Paper I and Paper II. However, this is not true in 
models with high-energy shocks since no extended ribbon 
is formed. These differences are apparent in Figure [5] 
which shows the final state of each model. Here we see 
that the mass in the high energy models is concentrated 
into a single object, moving outwards at a significant 
velocity from the dark matter halo. Conversely, for low 
E 55 values, the gas condenses to form a number of smaller 
star clusters that are well separated from each other. At 
the lowest energy level of E55 = 1, the shock strips some 
of the baryonic matter from the dark matter halo but 
leaves behind a small dense cluster at the center of the 
halo. Additionally, the metal content of such a weak 
shock is not sufficient to enrich the resulting cluster to 



the w 10 _2 Z© threshold. 

3.3. Effect of Minihalo Virialization Redshift 

Many of the basic properties of the halo, from the cen- 
tral density of the cloud to the radial density distribution, 
are strongly dependent on the minihalo virialization red- 
shift (see Eqs. HE]). In general, the higher the virializa- 
tion redshift, the more compact the minihalo, and thus 
the more resistant it is likely to be to external shocks. To 
explore the impact of this parameter, we conducted two 
simulations in which z c was taken to be 10 (PZC10) and 
15 (PZC15) as compared to our fiducial value of 8. The 
final states of each of these runs is shown in Figure El It 
is important to note that the higher redshift simulations 
act as a more robust prediction of our model for globular 
cluster formation since by z ~ 8 many of these minihalo 
may have been stripped of their baryons via ionization 
fronts during reionization. 

For all virialization redshifts, most of the mass from the 
halo is moved into a ribbon of material along the x-axis, 
but at higher virialization redshifts, the cloud is more 
compact (Eq. [T]) and the gas in the center of the halo 
cannot escape as easily as the surrounding gas. Also, 
as the virialization redshift increases, the physical size 
of the ribbon increases as does the total the entrained 
mass. Again, in all cases the clouds have been enriched 
to approximately 10 ~ 2 Z although the centers of the 
two oldest minihalos (those in runs PZC10 and PZC15) 
have slightly lower abundances. 

Figure [7] shows the evolved distribution of these mod- 
els. At lower virialization redshifts (z c = 8 and 10) the 
final distributions are very similar, and the number of 
clusters formed and the final mass of these clusters match 
well. At the highest virialization redshift, on the other 
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Fig. 3. — Comparison of clusters generated in models with vary- 
ing halo mass at a time of 200 Myrs after the end of the simulations. 
The top panel shows the logarithmic mass and the bottom panel 
shows the velocity of each cluster. The x-axis is the distance of 
each cluster from the center of their dark matter halos. The (red) 
unfilled squares show the PM03 model, the (blue) crosses show the 
OFID model, the (magenta) stars show the NFID/PM10 model, 
and the (black) filled squares show the PM30 model. The relative 
size of the points is proportional to the fraction of the total mini- 
halo gas mass contained in each cluster. Only clusters that have 
masses greater than 5% of the initial baryonic halo mass are shown. 
The total mass of clusters from each model is 4.1 xlO 4 , 4.0 xlO 5 , 
1.6 xlO 6 , and 4.8 xlO 6 M , for PM03, OFID, PM10, and PM30 
respectively. Except for the high mass-case, every cluster is found 
outside of its dark matter halo. 

hand, the distribution is quite different as one very large 
cluster is found nearly 20 kpc from the dark matter halo 
center. All of the other high redshift clusters are much 
smaller than those formed from minihalos with later viri- 
alization redshifts. This suggests that the formation of a 
large cluster (M & few xlO 5 M©) is robust over a large 
redshift range and may actually occur more easily for 
minihalos with earlier virialization redshifts. 

3.4. Effect of Shock Redshift 

Other than ^55, the most important parameter in de- 
termining the properties of the shock is its redshift, z s . 
As the shock redshift is increased the mass in the shock 
increases (Eq. [7)) and the metal abundance decreases 
(Eq. [TO]) , which in turn affect all aspects of the inter- 
action. To study the effect of a shock occurring at a 
higher redshift than our fiducial run, a model was run in 
which the outflow reached the minihalo at precisely the 
minihalo virialization redshift (z s = z c = 10). The final 
state of this run is contrasted with the fiducial case in 
Figure [8j 

The increased mass in the shock in run PZS10 leads 
to a much more stretched, elongated, and uniform mass 
distribution than in the NFID run. This elongated distri- 
bution makes it easier for the metals carried by the shock 
to be efficiently mixed into the collapsing cloud, and thus 



the metal abundance of the cloud in this run is higher 
than the fiducial model, even though the abundance of 
the shock is lower. On the other hand, the stretching in 
run PZS10 is not so severe as to suppress the formation 
of clusters, such as the one clearly visible at a distance 
of ^2.5 kpc from the center of the dark matter halo. 

Figure [9] compares the final stellar clusters generated 
by these two models. Even though the shock carries more 
mass in run PZS10 than in the fiducial case, the outcomes 
are fairly similar. Both have large clusters found outside 
the dark matter halo. Although fewer clusters are formed 
in run PZS10, the total mass in the clusters is slightly 
higher than in the fiducial model. 

3.5. Effect of Distance from Starburst Galaxy 

In Scannapieco et al (2004) and in Paper I and Pa- 
per II, the distance between the minihalo and the parent 

galaxy was taken to be R s = 1.5 M~ 1/6 (££ 55 ) (^f^)" 1 
kpc, a value that was chosen based on the observed 
metallicities of globular clusters. While this estimate 
gives a fiducial value of 3.6 (physical) kpc, in this study 
we consider the impact of varying this key parameter over 
a wide range of distances between R s = 2.1 kpc (PR21) 
and R 8 = 12 kpc (PR120). 

Figure [10] shows the final outcome from each of these 
models. When the minihalo is close to the source of 
the outflow, the cloud is quickly crushed into a dense 
cluster, similar to what occurs in runs with a larger shock 
energy. In fact, this happens so quickly that there is 
little time for the metals to be effectively mixed into this 
cluster and the gas ends with a lower metal abundance 
than our previous results. However, enough momentum 
is imparted to the gas to remove it from the dark matter 
halo. 

At intermediate distances (between 3.6 kpc and 6.6 
kpc) the outcome is familiar. The cloud is stretched into 
a ribbon of gas and is expelled from the dark matter halo. 
The more the distance increases, the more the cloud is 
stretched, and there is always enough time for metals to 
be mixed in the primordial gas to levels consistent with 
our previous runs. 

Finally, at the larget R s values, the cluster evolution is 
radically different than in the other models. The cloud 
is stretched as before but the resulting cluster does not 
leave the halo. The abundance of metals in the shock at 
this distance is already low, much lower than even the 
final cluster metallicities in many of our other runs, and 
while some metals are mixed into the cluster, it is very 
deficient compared to our other models. 

Figure [IT] shows the evolved state for each of these 
models. The panels are the same as Figure [5] The 
closest halo to the galaxy is transformed into a single 
cluster with a total mass of 7.0 x 1O 5 M . The interme- 
diate distance halos form large clusters with masses of 
w 3.0 x 1O 5 M and several other smaller clusters with 
velocities and positions consistent of being free from the 
dark matter halo. Finally, the farthest halo forms two 
dense clusters, however the largest cluster is found at 
the center of the dark matter halo, and only the smaller 
cluster is free of the halo potential at 5 kpc from its cen- 
ter. 

This suggests that there is a preferred distance from 
the outflow at which enriched clusters form most em- 
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Fig. 4. — Comparison of final outputs from runs with varying outflow energies. From top to bottom each row represents runs with .E55 
= 1, 5, 10, 20, and 30, respectively. As in Figure [2] the first column shows logarithmic density contours from 10 -26 to 10 -21 g cm -3 
and the second column shows logarithmic contours of metallicity from 10 -4 to 10 -0 ' 3 Zq, but now the third column shows logarithmic 
temperature contours between 10 and 10 8 ' 6 K. The length of the ribbon is correlated with the initial shock energy: the smaller the energy, 
the longer the ribbon. 



ciently. Too close to the outflow and the minihalo is 
crushed before it is enriched and too far from the outflow 
and the minihalo in neither enriched nor ejected from its 
dark matter halo. Therefore it seems that between « 3 
and 7 (physical) kpc is a preferable distance from a typ- 
ical starburst for the formation of compact stellar clus- 
ters. The observed distribution of halo globular cluster 
positions shows a drop-off beyond a galactocentric dis- 
tance of 40 kpc (Dauphole et al 1996). After evolving 
these clusters, the typical distance from the galaxy is be- 
tween 15-30 kpc, which agrees nicely with the observed 
distances. 



3.6. Effect of Halo Spin 

Next we explored the effect of net rotation on the evo- 
lution of the minihalo. In this case, the gas was given an 
initial velocity according to 



v x = -av c ^-cos((j)), 
H c 



v y av c 



x cos(</>) + z sin(</>) 
R c 



(11) 



v z = av c —sm(^>), 



where a is a constant, v c is the virial velocity of the 
cloud, R c is the virial radius, x and y are the positions 
within i? c , and <fi is the rotation angle. Here, a is set at 
0.05 and we vary <p between and 90 degrees so that the 
halo rotates around the z-axis and x-axis respectively. 
For this value of a we can calculate the spin parameter 



of our halos using the form from Bullock et al (2001b), 



A' 



J 



(12) 



where J is the angular momentum of the halo with mass 
M v i r contained in a sphere of radius R v i r and has a cir- 
cular velocity of V v i r . For the halos studied this gives a 
spin parameter of A' = 0.023, which is within 1 a of the 
mean value of X f = 0.035. 

Figure [12] compares the simulations at the final time. 
It is obvious that neither the spin direction nor the mag- 
nitude of the spin changes the final distribution. In all 
cases the cloud is stretched into a ribbon along the x- 
axis, the primordial gas is enriched to nearly constant 
value near 10 -2 Z , and the cluster is cooled to a few 
hundred degrees K. 

The reason for this insensitivity to X is a result of the 
shock itself. As the shock advances toward the halo it 
begins to develop vorticity, defined as uj = V x v, which 
evolves as ^ ~ ^Vp x Vp, where v is the velocity, 
p is the density, and p is the pressure. Physically, this 
baroclinic source term is a measure of the generation 
of vorticity due to the mismatch between the gradients 
of pressure and density (Glasner et al 1997). Since the 
shock is not aligned with the density gradient of the halo, 
this term is large throughout the simulation and vorticity 
begins to grow rapidly. 

Figure [T3] shows the z-component of vorticity as the 
shock heads toward the halo. By the time they meet the 
magnitude of the vorticity is much greater than the spin 
of the halo. In fact, the vorticity would be much greater 
than the spin of the halo even if we had chosen an a of 
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Fig. 5. — Comparison of clusters generated in runs with varying 
outflow energies. Panels are the same as Figure [3] except the y-axis 
of the top panel is given in linear units rather than logarithmic. The 
(green) unfilled triangles show the clusters formed from PE1, (red) 
unfilled squares show the clusters from PE5, the (blue) crosses show 
the clusters from NFID, the (magenta) stars show the clusters from 
PE20, and the (black) filled squares show the clusters from PE30. 
The total mass of the clusters formed is 9.0 xlO 6 , 2.3 xlO 6 , 1.6 
xlO 6 , 9.5 xlO 5 , and 8.2 xlO 5 M in PE1, PE5, NFID, PE20, and 
PE30 respectively. In general, the larger the outflow energy, the 
more the initial halo is crushed rather than forming a long stream 
of gas. This leads to the formation of fewer, larger clusters. 

1.0. 

3.7. Effect of Halo Concentration 

As the name suggests, the concentration parameter, c, 
describes in Eqj3] how the density is aggregated toward 
the center of the halo. Less concentrated halos have lower 
maximum central densities and shallower density profiles, 
while more concentrated halos have higher central gas 
densities and steeper density profiles at large radii. It is 
defined as 

c = R 2 oo/R s , (13) 

where R200 is the radius that corresponds to a density 
that is 200 times the critical density (Navarro et al 1997) 
and R s is the inner radius of the cloud. Typical values for 
NFW halos at this redshift are 4.8 (Madau et al 2001). 
However, Bullock et al (2001a) used high-resolution N- 
body simulations to study the dark matter halo den- 
sity distribution and found that for a given mass the 
range in concentration parameter can be fairly large, la 
A (log c) = 0.18. We therefore studied interactions with 
halos that cover this spread in concentration. 

Figure [TH compares models of the minihalo after alter- 
ing the concentration parameter. The top, middle, and 
bottom rows present the results of runs with c = 3.2, 
4.8, and 7.3 respectively. As the halo becomes more con- 
centrated, the metals have a harder time enriching the 
center of the halo, and becomes harder for the gas to be 
pushed out of the halo, which can be seen in the first 



column where the halo with the lowest concentration is 
farther away from the center of the halo. 

Figure fT5l shows the long term evolution of these mod- 
els. Generally there is little difference between these 
models, which all contain at least one large cluster with 
a mass few xlO 5 M . However, since there is more gas 
at near the edge of the less concentrated halos, the total 
mass of the ribbon is greater. Finally, in the most con- 
centrated halo there is a portion of the gas that fails to 
leave the dark matter halo. 

3.8. Summary 

Table [2] shows the final outcomes from all of our runs, 
including the total mass of clusters formed, and the basic 
properties of the largest cluster. Our findings can be 
summarized as follows: 

1. The mass of the initial minihalo is closely cor- 
related with the final mass of the compact clus- 
ters, such that the more massive the initial halo, 
the more massive the total distribution of clusters. 
However, as the halo gets more massive, not all the 
gas is removed from the halo, and a sizable cluster 
forms at the center of the dark matter potential. 
In all cases the gas is enriched to nearly identical 
metal abundances, and only the center of the most 
massive minihalo is slightly less enriched. 

2. The total energy in the outflow primarily affects the 
degree to which the final distribution is stretched. 
At low energies, the gas is shaped into a long rib- 
bon of material that is enriched at a lower level 
than the fiducial model. As the energy increases, 
the ribbon becomes smaller until the halo is com- 
pletely crushed, forming a single dense cluster with 
a similar metallicity as in the fiducial case. 

3. The collapse redshift of the minihalo directly affects 
the number and distribution of collapsed clusters. 
Halos that collapse early each produce a massive 
cluster that is found far from the dark matter halo. 
Halos that collapse later form a large number of 
lower mass clusters. Generally the clusters' metal 
abundances are independent of minihalo collapse 
redshift, but the earlier redshift runs yield clusters 
with slightly metal- deficient cores. 

4. If the redshift of the shock is increased, it entrains 
more mass, moves slower, and stretches the mini- 
halo gas more efficiently. At lower redshift s, the 
shock crushes the cloud more efficiently, and a more 
compact ribbon is formed. At higher redshifts, the 
cloud is enriched slightly more than the lower red- 
shift cloud, however both have an average abun- 
dance of w 10 -2 Z . 

5. The distance between the starburst galaxy and the 
minihalo has a dramatic impact on all properties 
of the final distribution of clusters. At small sepa- 
rations the minihalo is quickly crushed and is defi- 
cient in metals. At slightly larger distances, metals 
are found at levels similar to our fiducial model, 
and the cloud is stretched into a ribbon of material, 
which becomes longer in runs with larger distances. 
At very large distances the cloud is disrupted but 
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Fig. 6. — Comparison of final outputs from runs with different minihalo virialization redshifts. From left to right the columns show 
logarithmic contours of density, metallicity, and temperature with limits as in Figure |4] The top row shows the final state of run PZC15 
(z c = 15), the middle row shows the final state of run NFID (z c = 10), and the bottom row shows the final state of run PZC8 (z c = 8). 
In highest z c redshift run, multiple large clusters are formed while in the lower z c runs, only one primary cluster is formed. In all runs, 
the metal abundance remains roughly constant and near Z « 10 -2 Zq in the dense portions of the cloud, and the gas at the center of the 
dark matter halo is slightly deficient in metals. 



the outflow is too weak to move much of the gas out 
from the dark matter halo. Because the shock itself 
is metal-poor, the final cluster is also metal-poor. 

6. The initial spin of the halo has no discernible im- 
pact on the evolution of the minihalo. This is 
because the vorticity generated during the shock 
minihalo interaction dwarfs that of the minihalo by 
orders of magnitude. All cases studied were iden- 
tical to the case without spin. 

7. The concentration of the halo has a strong effect on 
the positions of the collapsed stellar clusters. The 
more concentrated the minihalo, the harder it is for 
the outflow to remove the gas from the dark mat- 
ter potential, and in the most concentrated case we 
studied, some halo gas remains in the dark matter 
halo. In all cases the gas is enriched to nearly iden- 
tical levels with only the most concentrated gas at 
the center of the halo being slightly metal deficient. 
However, dense clusters are also always formed out- 
side of the dark matter potential. 

4. OBSERVATIONAL SIGNATURES 

The suite of simulations described above shows that 
over an extremely wide range of parameters, the ulti- 
mate consequence of the interaction between a galactic 
outflow and a primordial minihalo is the rapid forma- 
tion of dense clusters containing up to a few 1O 6 M of 
stars. While such high-redshift clusters are not directly 
observable with current telescopes, their rapid bursts of 
star formation and consequent low mass-to-light ratios 



present a opportunity for study with the next generation 
of instruments. At the same time, their compact na- 
ture and formation in low-density environments makes 
it likely that many of them may have survived to the 
present, allowing for indirect connections with current 
stellar populations. Here we explore both of these con- 
nections. 

4.1. Direct Observations 

To calculate the direct observability of the clusters 
formed in our simulations, we first constructed an es- 
timate of the number of stars formed as a function and 
time and position. As there is no explicit prescription for 
star formation included in our simulations, we instead 
built up the star-formation history by post-processing 
our outputs, carrying out the following steps: 

1. First we calculated the total stellar mass in a series 
of 175 evenly-spaced bins, adding together the mass 
in cells that exceeded a density threshold of p t h = 
1.0 x 10 -23 g cm -3 . Note that as the mass collapses 
dramatically onto the x axis, its final density is 
limited by the resolution in the y and z directions. 
Thus the threshold is less than the density of the 
collapsed cloud in nature, and it was determined by 
examining the late stages of our fiducial model and 
finding the apparent edge of the ribbon of material. 
In this way the collapsed mass was computed for 
every simulation output, such that we computed 
the masses and position of the stars as a function 
of time throughout the simulation. It should be 
noted that the masses quoted below are assume 
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TABLE 2 

Summary of model outcomes. M^aio is the initial gas mass of the minihalo in units of 10 6 M@, M c i usters is the total mass 

FOUND IN CLUSTERS IN UNITS OF 10 6 Mq, Mi arge IS THE MASS OF THE LARGEST CLUSTER FORMED IN UNITS OF 10 6 Mq, Z c i uster IS THE 
METALLICITY OF THE LARGEST CLUSTER, V c i uste r IS THE VELOCITY OF THE CLUSTER AFTER 200 MYRS OF EVOLUTION IN UNITS OF KM S _1 , 
AND D c i uster IS THE DISTANCE OF THE LARGEST CLUSTER FROM THE CENTER OF ITS DARK MATTER HALO IN UNITS OF (PHYSICAL) KPC. 



Name 


M ha l 


M c l us ters 




Zciuster 


V c l us ter 


-^cluster 


Notes 




10 6 


10 6 Mq 


10 6 Mq 


z 


l —1 

km s 


kpc 




OFID 


0.45 


0.40 


0.07 


10 z - u 


45 


9.8 




PM10/NFID 


1.5 


1.6 


0.38 


10 2U 


59 


13 




PM03 


0.05 


0.01 


0.01 


10 2U 


76 


16 




PM30 


4.5 


4.8 


1.9 


1Q -2.0 


63 


15 




PE1 


1.5 


0.9 


2.8 


LO-3.0 


17 


5.0 


Data for cluster outside of DM halo 


PE5 


1.5 


2.3 


0.37 


LO-2.5 


72 


17 




PE20 


1.5 


0.9 


0.66 


1Q -2.2 


41 


9.6 




PE30 


1.5 


0.8 


0.82 


10 -2.0 


53 


12 




PZC8 


1.5 


1.2 


0.32 


1Q -2.0 


60 


13 




PZC15 


1.5 


2.2 


0.96 


LO-2.0 


130 


29 




PZS10 


1.5 


2.3 


0.37 


1Q -1.5 


72 


18 




PR21 


1.5 


0.7 


0.64 


1Q -3.0 


41 


8.9 




PR66 


1.5 


2.1 


0.48 


1Q -2.0 


45 


13 




PR120 


1.5 


0.75 


0.19 


1Q -3.0 


14 


5.0 


Data for cluster outside of DM halo 


PC32 


1.5 


1.8 


0.54 


LO-2.0 


62 


14 




PC73 


1.5 


1.5 


0.28 


10 -2.0 


58 


13 


Data for cluster outside of DM halo 
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Fig. 7. — Comparison of clusters generated in runs with different 
virialization redshifts. Panels are the same as Figure El The (red) 
unfilled squares show the clusters from PZC8, the (blue) crosses 
show the clusters from NFID, and the (magenta) stars show the 
clusters from PZC15. The total mass of clusters is 1.2 xlO 6 , 1.6 
xlO 6 , and 2.3 xlO 6 M© for runs PZC8, NFID, and PZC15 respec- 
tively. At lower redshifts the sizes and distribution of clusters is 
fairly similar, with at least one cluster with mass greater than 2.0 
xlO^M© while at high redshift only one large cluster is formed. 
In all cases the largest clusters are found far from their respective 
dark matter halos. 

that the star formation efficiency is 100 percent, 
and therefore the quoted stellar mass and stellar 
flux are an upper limits. 

2. Next, we used the positions and velocities to corre- 



late the stellar distribution at each output with the 
distribution at the previous output. Starting with 
with the final output, and working back through 
the files in this way, we calculated the mass of new 
stars formed as a function of time and determined 
the positions of these stars at the final output time. 

3. Finally, we calculated the stellar fraction as a func- 
tion of age for each of the mass bins at the final 
output time. 

Using this information, we could then estimate the 
emitted flux for an arbitrary set of narrow or wide band 
filters. At redshifts z > 8 the James Webb Space Tele- 
scope (JWST) will be the best telescope to see these ob- 
jects in wide band filters. From the ground, the best 
method of detection is to search for redshifted Lyman 
alpha (Lya) emission using near-infrared capabilities of 
next generation, 30-40 meter class telescopes such as 
the Giant Magellan Telescope (GMT), the Thirty Me- 
ter Telescope (TMT), and the European Extremely Large 
Telescope (E-ELT). 

To determine the broad-band fluxes of our stellar 
clusters, we used the population-synthesis code bc03 
(Bruzual & Chariot 2003) to compute the luminosity 
per solar mass of a stellar population as a function of 
frequency and age, and convolved this with the stellar 
history of each bin, and related this to the flux at an 
observed frequency vq as 



(1 + *) 



dt 



dM* 



dM* 
~dt 



(*), (14) 



in units of ergs cm -1 Hz -1 s _1 and where z is the shock 
redshift, d\[z) is the luminosity distance, vq/(1-\-z) is the 
rest-frame frequency, c and dLv j dM*(y§j (\ + z), t) is the 
luminosity per frequency per solar mass of a population 
of stars with an age t, and dM*/dt(t) is the star formation 
history in a given bin. 
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Fig. 8. — Comparison of final outputs from runs NFID (top) and PZS10 (bottom) which have outflow redshifts of z s = 8 and z s = 10, 
respectively. Panels are the same as Figure |4] The density is much more stretched and uniform than the fiducial case. While the metals 
from the shock have been mixed into the primordial gas, it is significantly more enriched than in the fiducial model. Both models reach 
the same final temperature. 

photons and is taken to be 0.2, Q(H) is the hydrogen ion- 
izing photon rate in units of number per second per solar 
mass and is taken from the stellar population synthesis 
code STARBURST99 (Leitherer et al 1999). Finally, in 
both the broad and narrow band cases we normalized 
the flux by the spacing of each bin to get the flux per 
(physical) kpc. 

Figure [16] shows the resulting fluxes for three different 
JWST bands (F115W, F150W, and F200W) and the Ly- 
man alpha flux from a selection of simulations spanning 
a wide range of parameter space. Again we note that 
the fluxes in this figure are computed assuming that all 
the gas collapsed above p t h = 1.0 x 10 -23 g cm -3 is con- 
verted into stars, and so these fluxes should be scaled 
by an unknown star formation effciency factor which is 
< 1 . Likewise the angular given in this figure assume that 
that stellar clusters are being viewed edge on, such that 
angular separation is maximal. Table [3] summarizes the 
observable properties of the stellar clusters generated in 
all of our models including the total fluxes in the JWST 
bands and Lya and the physical and maximal angular 
scales. Note that the physical and angular scales given 
are for the extended emission expected at the end of each 
simulation and not from the final globular cluster. 

From the values in this figure and table we see that, 
unfortunately even the wide band filters of JWST are not 
expected to have the required sensitivity to detect these 
objects. Observations with JWST wide band filters will 
have typical sensitivities of 10-20 nJy for a 10-a detec- 
tion at 10,000 seconds of integration time (Stiavelli et 
al 2008), which is roughly an order of magnitude higher 
than the typical fluxes of « 1 nJy (see Figure [T6|) ex- 
pected from our objects. 

On the other hand, the narrow band Lya fluxes are 
over an order of magnitude brighter than those expected 
to be obtained with the next generation of ground- 
based telescopes. The proposed Near Infrared Multi- 
object Spectrograph (NIRMOS) on the GMT for exam- 
ple, will be able to detect such sources down to a flux 
limit of 1.0 xlO -20 ergs/s/cm -2 , given 25 hours of in- 
tegration time (McCarthy 2008; GMT Science Case). 
Similarly, the InfraRed Imaging Spectrometer (IRIS) on 



o - 



J I I I I I I I I I I I I I I I I I I I L 



5 10 15 

Position [kpc] 



20 



Fig. 9. — Comparison of clusters generated by runs with varying 
outflow redshifts. The panels are the same as in Figure [5] The 
(red) unfilled squares show the clusters from NFID and the (blue) 
crosses show the clusters from PZS10. The total mass of clusters for 
these models is 1.6 xlO 6 , and 2.3 xlO 6 M for NFID and PZC10 
respectively. Although the shock lasts longer in the PZS10 case, 
there is little difference between these two cases. The PZS10 model 
has fewer, but slightly larger clusters than in the fiducial case and 
they are found outside of its dark matter halo. The positions and 
velocities are comparable between these two models. 

Similarly, we estimate the Lyman alpha flux from our 
simulations as 



1 



4*4} (z) 



I 



dL a dM* 

dt ^-r{t)——{t), 



dM* 



dt 



(15) 



where dL a /dM*(t) = c L (l - f esc )Q(H)M*, is the the Ly- 
man alpha luminosity per solar mass of a population of 
stars with an age t (e.g. Scannapieco et al 2003), with cl 
= 1.04 xlO -11 ergs, f esc is the escape fraction of ionizing 
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Fig. 10. — Comparison of final outputs from runs with varying distances between the minihalo and the galaxy. The first column shows 
logarithmic density contours between 10 — 26 and 10 — 21 g cm -3 , the second column shows logarithmic metallicity contours between 10 — 4 
and 10 -0-3 Zq, and the third column shows logarithmic temperature contours between 10 and 10 8 6 K. The top row shows the result of 
placing the minihalo at a distance of 2.1 kpc (PR21), the second is the fiducial distance of 3.6 kpc (NFID), the third row is a model with 
a distance of 6.6 kpc (PR66), and the last row shows the farthest case with a distance of 12.0 kpc (PR120). Beyond the apparent physical 
differences between each model, the metal abundance is lower than our fiducial model if the halo starts too close or too far from the outflow. 
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Fig. 11. — Comparison of clusters generated by runs with varying 
distances between the minihalo and the galaxy. The (red) unfilled 
squares show the clusters formed in PR21, the (blue) crosses show 
the clusters formed in NFID, the (magenta) stars show the clusters 
formed in PR66, and the (black) filled squares show the clusters 
formed in PR120. The panels are the same as in FigureO The total 
mass in each model is 7.0x 10 5 , 1.6x 10 6 , 2.1 x 10 6 , and 7.5 x 1O 5 M 
for PR21, NFID, PR66, and PR120 respectively. While all models 
create at least one dense cluster far from the center of the halo, 
at the largest distances from the starbursting galaxy, the outflow 
is not strong enough to remove all of the gas from the halo. This 
leaves a cluster at the center of the dark matter halo. 



the TMTwill detect Lya sources at z = 7.7 with fluxes 
of 1.0 xlO -18 ergs/s/cm -2 with signal-to- noise (S/N) of 
15 in only 1 hour of observation time (Wright & Bar- 
ton 2009; TMT Instrumentation and Performance Hand- 
book). Finally, OPTIMOS-EVE (Optical- Near- Infrared 
Multi-object Spectrograph) for the E-ELT will detect 
sources with fluxes of 10 -19 erg/s/cm -2 with S/N of 8 in 
40 hours of integration time (Hammer et al 2010). This 
means that the majority of the models studied here are 
more than bright enough to be detected. 

Beyond being bright in Lya, the unusual, elongated 
morphology of the stellar distributions formed in outflow- 
minihalo interactions makes them ideal for study with 
next generation of ground-based instruments. With ex- 
pected angular resolutions of 0.1 — 0.3 arcsec (McCarthy 
2008; Wright & Barton 2009; Hammer et al 2010), next 
generation narrow band images will not only be able to 
detect the presence of the stars, but show that their dis- 
tribution is highly-elongated in the direction of the im- 
pinging outflow, as illustrated in Figure[17l Furthermore, 
as the starburst galaxies triggering star-formation should 
have typical solar masses J> 1O 8 M , these will be easily 
detectable with broad-band J WST measurements. Thus, 
the detection of a group of elongated Lya emitters which 
whose axes point directly at a larger broad-band detected 
starbursting galaxy will provide a unique signature that 
unambiguously points to the formation of compact stel- 
lar clusters by high-redshift galaxy outflows. 

4.2. Connection with Halo Globular Clusters 

Equally exciting as the prospect of direct detection 
is the connection between outflow-induced star forma- 
tion at high-redshift and stellar clusters observed around 
nearly galaxies. Over a large range of parameter space, 
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Fig. 12. — Final profiles of the spinning versus non-spinning halos. The top row shows the case in which the spin is about the z-axis 
(PSPZ). The middle row shows the non-rotating case (NFID) and the bottom row shows the case where the spin is about the x-axis 
(PSPN). The first column shows the logarithmic density contours, the second column shows the logarithmic metallicity contours, and the 
third column shows the logarithmic temperature contours. Panel limits are the same as in Figure [2] There is very little difference between 
each of these runs, which suggests that the initial spin of the minihalo does not contribute significantly to its final evolution. 
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Fig. 13. — Snapshots of the ^-component of vorticity at different 
times before the shock impacts the halo in run PSPZ. The top left 
panel shows the initial vorticity of the cloud. The top right and 
both bottom panels show the evolution of vorticity as the shock 
nears the halo. The scale is in units of km/s/kpc. As the shock 
nears the minihalo the vorticity increases and by the time it reaches 
the halo, it is much larger than the spin of the halo. 



the interactions studied in our simulations result in com- 
pact, tighty-bound stellar clusters that are found in low- 
density environments in which cluster disruption pro- 
cesses are minimal. Thus it is likely that a substantial 
fraction of the clusters generated by outflow-minihalo in- 
teractions may persist even to the present day. As men- 
tioned above and in our previous investigations (Paper I 
and Paper II), a natural low-redshift counterpart of these 
compact, high-redshift clusters may be the population of 



TABLE 3 

Summary of simulated model fluxes. F115W is the total 

FLUX FROM THE F115W JWST BAND IN UNITS OF NjY, LY a IS 
THE TOTAL LYMAN OL FLUX IN UNITS OF ERGS/s/CM 2 , AND THE 
THIRD AND FOURTH COLUMNS ARE THE MAXIMUM ANGULAR AND 
PHYSICAL SPATIAL SCALES RESPECTIVELY. 



Model 


F115W 




arcsec (") 


kpc 


OFID 


0.25 


9.2e-19 


0.27 


1.3 


NFID 


1.02 


2.8e-18 


0.28 


1.4 


PM03 


0.02 


2.3e-20 


0.14 


0.7 


PM30 


2.7 


3.0e-18 


0.55 


2.7 


PE5 


1.3 


1.2e-18 


0.65 


3.2 


PE20 


0.54 


6.1e-19 


0.16 


0.76 


PE30 


0.40 


4.6e-19 


0.08 


0.39 


PZC8 


0.53 


5.4e-19 


0.24 


1.2 


PZC15 


2.1 


2.3e-18 


0.46 


2.3 


PZS10 


0.89 


8.9e-19 


0.68 


2.9 


PR21 


0.35 


4.0e-19 


0.05 


0.25 


PR66 


0.47 


2.3e-19 


0.54 


2.7 


PC32 


0.99 


8.5e-19 


0.27 


1.3 


PC72 


0.68 


7.2e-19 


0.30 


1.5 



halo globular clusters. 

How globular clusters are formed has been an open 
question since their discovery. The typical age of a glob- 
ular cluster is between 10-13 Gyrs (Krauss & Chaboyer 
2003), which demands a cosmological origin for these sys- 
tems. Observations of galactic globular clusters have 
found a bimodal color distribution which suggests two 
subpopulations (Zepf & Ashman 1993). The relatively 
metal rich population with \^r] ~ —0.5 are associ- 
ated with the thick disk or bulge of galaxies and are 
thought to have formed as a consequence of galaxy in- 
teractions during the build-up of their parent galaxy 
(Shapiro et al 2010; and see Brodie & Strader 2006 for 
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Fig. 14. — Comparison of final outputs from runs with different concentration parameters. The first row shows the c = 3.2 case (PC32), 
the second row shows the fiducial c = 4.8 case (NFID), and the third row shows the c = 7.3 case (PC72). Columns are the same as Figure 
[2] The more concentrated the halo, the more stretched it becomes and the less it is enriched. 
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Fig. 15. — Comparison of clusters generated in runs with differ- 
ent halo concentrations. Panels are the same as Figure [5] The 
(red) unfilled squares show clusters from the least concentrated 
halo (PC32), the (blue) crosses show the clusters from the fiducial 
cluster (NFID), and the (black) stars show the clusters from the 
most concentrated halo (PC73). The total mass of clusters in each 
model is 1.8 xlO 6 , 1.6 xlO 6 , and 1.5 xlO 6 M for PC32, NFID, 
and PC73 respectively. The outflow has a hard time removing the 
gas from the most highly concentrated halo, and in this run a clus- 
ter is formed within the dark matter potential well. On the other 
hand, all the gas is ejected from the least concentrated halo, which 
forms the most massive cluster of any of the three runs. 

a review). On the other hand, the metal-poor popula- 
tion with [^J « — 1.6 are associated with galaxy halos 
and are thought to form in the early universe through an 
unknown mechanism. 



Formation scenarios for such clusters are typically split 
into two groups: 'pre-enrichment' and 'self-enrichment' 
schemes. In the pre-enrichment picture, the primordial 
gas is homogeneously enriched via supernovae in such 
a way that does not disrupt the cloud (e.g. Beasely et 
al 2003; Bromm & Clarke 2002; Elmegreen & Efremov 
1997; Marcolini et al 2009). However, little is known 
about this previous generation of stars, how it was able 
to enrich the gas so quickly, and why it played only a sec- 
ondary role in the formation of globular clusters. In the 
self-enrichment picture the primordial gas was enriched 
by a supernova contained within the halo hosting the 
forming globular cluster (e.g. Boley et al 2009; Brown et 
al 1995; Recchi et al 2005; Smith 2010). This has the 
problem of not being able to mix the metals into the 
cluster rapidly enough, as well as usually leads to the 
the supernova unbinding the cluster (Peng & Weisheit 
1991; Whalen et al 2008b). Finally, Cen (2001) set aside 
the question of enrichment completely and used a sim- 
ple model to suggest that ionization fronts can act as an 
external force to collapse these clouds. However, many 
others (e.g. Haiman et al 2001; Illiev et al 2005; Shapiro 
et al 2004) showed that instead of forming a dense clus- 
ter, the ionization front completely boils away the mini- 
halo. As shown in our simulations, interactions between 
galaxy outflows and minihalos share the best aspects of 
all these scenarios: bringing in metals only moments be- 
fore vigorous star formation commences, depositing the 
metals without unbinding the proto-globular cluster, and 
triggering collapse without evaporating the cloud. 

Observationally, there are three important properties 
of halo globular clusters that make this connection a 
promising one. First, globular cluster masses are dis- 
tributed as a Gaussian with a mean value of 1O 5 M and 
a dispersion of 0.5 dex (e.g. Armandroff 1989). Here, the 
lower mass cut-off in this population is likely to be due 



16 



0.0 0.5 1.0 1.5 



Position [kpc] 

0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 



o 

Q. 



2.0 
1.5 
1.0 
0.5 



o 

Q. 



O 



2.0e-18 
1.5e-18 
1.0e-18 
5.0e-19 

1.00 
0.80 
0.60 
0.40 
0.20 



o 

Q. 



O 
Q. 



O 

w 



8.0e-19 
6.0e-19 
4.0e-19 
2.0e-19 




OFID 





R s =6.6 




0!4 0:9 1:3 17 2M 2.6 
Z c =10/ 




j i i_ 



0.0 0.1 0.2 0.3 



0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 
Position ["] 



Fig. 16. — Simulated fluxes from selected models, scaled to 100% star formation efficiency. The top two rows show the simulated fluxes 
for runs NFID, PE5, and PE20 (from left to right), while the bottom two rows shows the simulated fluxes for runs OFID, PR66, and 
PZS10. The first and third rows show the expected fluxes in JWST wide band NIRCam filters. The (red) solid lines are fluxes in the 
F115W band, the (dashed) green lines are in the F150W band, and the (dotted) blue lines are in the F200W band, all measured in njy per 
kpc. The second and fourth rows show the expected Lyo; intensities and are measured in ergs/s/cm 2 /kpc. The top x-axis is the physical 
spatial scale measured in kpc of each model while the bottom x-axis is the angular scale measured in arcseconds. The y-axis is identical 
across a given row. 



to a variety of destruction processes including mechani- 
cal evaporation {e.g. Spitzer & Thuan 1972) and shock- 
ing as the cluster moves through the disk of the galaxy 
(e.g. Ostriker et al 1972), but the high-mass cut-off ap- 
pears to be a property of the initial population. Except 
for the models with the smallest minihalo masses, all 
of our models produce at least one dense cluster with a 
mass between 10 5 and a few times 10 6 M . Furthermore 
minihalos have an intrinsic maximum mass of « 10 7 M 
which corresponds to the T w 10 4 K limit where atomic 
hydrogen/helium cooling becomes inefficient, and thus 
the maximum sizes of compact clusters in our study are 
likely to place a rough upper bound on the masses of 
stellar clusters that can be formed by this mechanism in 
nature. 

A second important property connecting our high- 
redshift clusters with the present-day population of halo 
globuar clusters is the abundances of stars both within 
a given cluster and between different halo globular clus- 
ters. The metallicity distribution between clusters is well 
defined by a Gaussian distribution with a mean value 
[if] ~ —1-6 with a dispersion of 0.3 dex (Zinn 1985; 



Ashman & Bird 1993). Most individual clusters have a 
dispersion of less than 0.1 dex (see Suntzeff 1993 and ref- 
erences within) , although it is worth noting that some of 
the clusters that show larger scatter may be due to sub- 
sequent star formation and enrichment due to evolved 
stars (e.g. Piotto et al 2007; D'Ercole et al 2008; Bekki 
2011). 

In our simulations, the metals from the incoming out- 
flows are well mixed into the primordial gas from the 
minihalo through the turbulent processes that are inher- 
ent to this interaction. Over a wide range of parameters, 
the metal-free gas is enriched to nearly constant value of 
Z w 10 -2 Z , which approximately matches the observa- 
tions. Only a small subset of models are found where the 
metal abundance is below this value. In cases where the 
halo is very close to the galaxy, the minihalo is crushed 
before an appreciable abundance of metals is transported 
into the primordial gas. If the halo is too far away from 
the galaxy then there is ample time for metals to move 
into the cloud, but the shock is too deficient to enrich 
to the fiducial level. In most models there is some slight 
difference in the metal abundance especially for the gas 
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Fig. 17. — Simulated narrow-band images from selected models. The top two rows show simulated images for runs NFID, PE5, and 
PE20 (from left to right), while the bottom two rows shows the simulated images for runs OFID, PR66, and PZS10. The first and third 
rows show forming clusters as observed edge-on with 0.25 arcsecond resolution, and the second and fourth rows show clusters as observed 
edge-on with 0.1 arcsecond resolution. In each panel the x and y axes are in units of arcsec, and the contours are labeled according to flux 
per unit area relative to the maximum flux per unit area in the image. 



originally found near the center of the dark matter halo. 
Otherwise, it seems that metal enrichment in these sit- 
uations is fairly uniform and robust around a value of 

«10" 2 z . 

Thirdly, the observation that globular clusters do not 
reside within dark matter halos provides a strong con- 
straint on their formation. Such observations show that 
tidal forces are actively stripping more stars from glob- 
ular clusters (Irwin & Hatzidimitriou 1993; Grillmair et 
al 1995) than would be expected if housed within dark 
matter halos (Moore 1996; Conroy et al 2010). This too 
is a robust prediction of our model, and in only in a small 
subset is there any gas left in the dark matter halo. In 
fact, only in runs with the largest halo mass, the largest 
concentration, and at the largest separations between 
minihalo and the starburst galaxy was gas retained by 
the minihalo, and in all these cases at least one other 



cluster was formed that was unbound from the halo. 
Thus outflow-minihalo interactions are a mechanism that 
primarily, but not exclusively produces dark-matter free 
clusters. Presumably these rare clusters would still be 
found with their dark matter halos if they have not been 
stripped away by some other means or buried within the 
center of large, low-redshift galaxies. 

Globular cluster formation has also been studied using 
simulations of the heirarchical buildup of a Milky Way 
sized galaxy. Kravtsov et al (2005) carried out one such 
a simulation and found that proto-globular clusters are 
produced in giant molecular clouds within the disk of 
the galaxy. While this model reproduced many of the 
properties expected of halo globular clusters, it relied 
on subsequent violent mergers to move these clusters to 
the galaxy halo. Muratov & Gnedin (2010) looked at a 
similar mechanism that reproduces the observed metal- 
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licity distributions found in globular clusters. GrifTen 
et al (2010) examined a halo from the Aquarius simula- 
tion for sites of globular cluster formation, and adopted 
a simple model in which the cluster sites are determined 
solely on their temperature. Cluster formation is ended 
with when the host galaxy is completely reionized. They 
were able to reproduce the expected number of present 
day clusters, their positions, and formation ages. How- 
ever, since their simulation uses only dark matter parti- 
cles, no mention is made of how the baryonic matter is 
removed from the dark matter halo or how the metal- 
licity of the resulting cluster arises. On the other hand, 
the minihalo-galaxy interactions studied here are both 
a natural consequence of heirarchical galaxy formation, 
and they reproduce the masses, metallicities, and dark 
matter content of halo globular clusters directly in our 
simulations. 

5. CONCLUSIONS 

The early universe was permeated by primordial mini- 
halos that provided the building blocks for larger, later- 
forming structures. However, since they were not mas- 
sive enough to form stars on their own, these objects 
remained passive until acted on by an outside influ- 
ence. This means the first galaxies were formed in some- 
what larger dark- matter halos, which cooled atomically, 
formed stars and supernovae, and funned a fraction of 
the resulting energy into massive galaxy-sized outflows. 
For minihalos in orbit around these early galaxies, such 
outflows may have triggered a radical transformation. 

In our previous papers, we studied the role of pri- 
mordial chemistry and cooling (Paper I) and turbulent 
mixing and metal-line cooling (Paper II) on outflow- 
minihalos interactions with a fiducial set of parameters. 
In each of these simulations, the baryonic matter was 
expelled from the minihalo potential and formed dense, 
cold clusters with nearly constant metal abundances. 
Noting that these interactions could be important in un- 
derstanding the origin of halo globular clusters, we care- 
fully studied each aspect of this interaction and showed 
that this result was independent of uncertainties in the 
chemistry and turbulence models and the level of the 
high-redshift background of dissociating photons. Fur- 
thermore, varying the maximum resolution in our simula- 
tions, we showed that medium-resolution (256 x 128 x 128 
effective) simulations were able to faithfully reproduce 
the outcome of these collisions. 

Here we complete this picture and perform a large, 
medium-resolution parameter study, quantifying the im- 
pact of minihalo mass, minihalo formation redshift, out- 
flow energy, outflow redshift, distance, minihalo concen- 
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tration, and spin. For a wide range of parameters, the 
results are extremely similar. The baryonic matter is 
expelled from the dark matter halo and formed into at 
least one dense, cold cluster that is homogeneously en- 
riched with metals. In fact only under extreme circum- 
stances, such as a large separation between the halo and 
the galaxy, very low energy outflows, or very high mini- 
halo concentration, is gas retained by the minihalo, and 
even in these cases at least one other compact, unbound 
cluster is formed. 

Furthermore, our parameter study strengthens the idea 
that the longest-lived stars formed by this processes will 
be observable today as members of halo globular clus- 
ters. Like the clusters in our simulations, such globular 
clusters are observed over a substantial mass range, and 
their upper mass limit can be directly associated with the 
maximum minihalo mass, above which atomic cooling be- 
comes efficient. Over a wide range of energies, redshift s, 
and distances, outflows are able to accomplish three im- 
portant jobs necessary to form realistic halo globular 
clusters: imparting the momentum required to move the 
pristine gas from the dark matter halo, starting the non- 
equilibrium chemistry and cooling required for collapse, 
and providing a source of metals. The turbulence that 
follows then mixes these metals into the primordial gas 
nearly homogeneously. 

While the direct detection of outflow-minihalo interac- 
tions is beyond current capabilities, it will be well within 
the reach of the telescopes currently being planned. Post- 
processing our simulations we show that the outflow- 
driving galaxies are likely to be detectable in broad-band 
JWST images, but the clusters themselves are likely to be 
just beyond their expected detection limits. On the other 
hand, narrow-bound imaging of redshifted Lyman alpha 
emission from these forming clusters will be well within 
the capabilities of large ground-based telescopes like the 
GMT, TMT, and E-ELT. Such Lya emitters will appear 
as bright, extended in a single direction, and pointed di- 
rectly at larger broad-band detectable starbursts. This 
unique signature makes them perfect targets for the next 
generation of telescopes and an exciting observational 
probe of an extraordinary mode of high-redshift star for- 
mation. 
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